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Abstract: We present new parameterizations of vector and axial nucleon form factors. We maintain an 
excellent descriptions of the form factors at low momentum transfers, where the spatial structure of the 
nucleon is important, and use the Nachtman scaling variable £ to relate elastic and inelastic form fac- 
tors and impose quark-hadron duality constraints at high momentum transfers where the quark structure 
dominates. We use the new vector form factors to re-extract updated values of the axial form factor from 
neutrino experiments on deuterium. We obtain an updated world average value from neutrino and pion 
electroproduction experiments of Ma = 1.0144 ± 0.0136 GeV/c 2 . Our parameterizations are useful in 
modeling atmospheric neutrino interactions (e.g. for neutrino oscillations experiments). 
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At low Q 2 , a reasonable description of the pro- 
ton and neutron elastic form factors is given by the 
dipole approximation. 
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where C V ' A = (l,g A ), 9a = -1-267, My = 0.71 
(GeV/c 2 ) 2 , and M A = 1.015 GeV/c 2 . 

Here we present parameterizations that simultane- 
ously satisfy constraints at low Q 2 where the spa- 
tial structure of the nucleon is important, as well 
as at high Q 2 where the quark structure is impor- 
tant. Our new quark-hadron duality based param- 
eterization ("BBBA07") includes: (1) Improved 
functional form that uses Nachtman scaling vari- 
able £ to relate elastic and inelastic vector and ax- 
ial nucleon form factors; (2) Yield the same val- 
ues as fits by Arlington and Sick [3] for Q 2 < 
0.64(GeV '/ 'c) 2 , while satisfying quark-hadron du- 
ality constraints at high-Q 2 . 

We use our new "BBBA07"vector form factors to 
re-extract updated values of the axial form factor 
from a re-analysis of previous neutrino scattering 
data on deuterium and present a new parameteri- 
zation for the axial form factor within the frame- 
work of quark-hadron duality. For vector form fac- 
tors our fit functions are A/v(£) (i.e. Ae p (£), 



A-Mp{£,), A En (£), A Mn (0) multiplied by an up- 
dated Kelly [3] type parameterization of one of the 
proton form factors. 
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where clq = 1, m = 1, and r = Q 2 /4M^, 
(Af/v is proton, neutron, or average nucleon mass 
for proton, neutron, and axial form factors, respec- 
tively). The datasets used by Kelly to fit Ge p and 
Gmp/Hp (Mp = 2.7928, p n = —1.913) are de- 
scribed in [3], Our parameterization employs the 
as-published Kelly parameterization to G^p lly and 

an updated set of parameters for G^ e p V (Q 2 ) that 
includes the recent BLAST[2] results. The param- 
eters used for G Kelly (Q 2 ) are given in ref. [1]. 
Each Pj is a LaGrange polynomial in the Nacht- 



man variable, £ 



The £j are equidis- 



(l + Vl+l/r) 

tant "nodes" on an interval [0, 1] and pj are the fit 
parameters that have the property Apf(£j) = Pj. 
The seven pj parameters are at £j=0, 1/6, 1/3, 1/2, 
2/3,5/6, and 1.0. 
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Figure 1: Ratios of Ge p (a), Gm p /h p (b), Geu (c) and Gm„/ 'fi n (d) to G^. The short-dashed line in 
each plot is the old Kelly parameterizations (old Galster for Geu)- The solid line is our new BBBA07 
parameterization for ^ = 0.0, and the long-dashed line is BBBA07 for ^ = 0.2. 



In the fitting procedure described below, the pa- 
rameters of Am (£) are constrained to give the same 
vector form factors as the recent low Q 2 fit of 
Arrington and Sick [3] for Q 2 < 0M(GeV/c) 2 
(as that analysis includes coulombs corrections 
which modify G£p,two photon exchange correc- 
tions which modify Gm p and Gmu)- Since the 
published form factor data do not have these cor- 
rections, this constraint is implemented by in- 
cluding additional "fake" data points for Q 2 < p 2 — e 2 qi (£)) 



values A(C)=Px at £ 1= (Q 2 = 0) for G Mp , G Ep , 
Geui Gun ar e set to to 1.0. The value A(^)=p- 7 at 
£j = l (Q 2 — * °°) for Gm p and Ge p is set to 1.0. 
The value A(£)=pj at £,=1 for Gmu and Geu are 
fixed by constraints from quark-hadron duality. 
Quark-hadron duality implies that the ratio of neu- 
tron and proton magnetic form factors should be 
the same as the ratio of the corresponding inelas- 
tic structure functions in the f=l limit. (Here 



0.64(GeV/c) . Our fits to the form factors are 
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where we use our updated parameters in the Kelly 
parameterization of Gm p - For Geu the parameters 
a=1.7 and b=3.3 are the same as in the Galster[3] 
parametrization and ensure that dGEn/dQ 2 at for 
Q 2 = is in agreement with measurements. The 



We ran fits with two different values of - at the 
£=1 limit: ^ =0 and 0.2 (corresponding to = 
0.25 and 0.4286). The fit utilizing | = is Af In , 
and the fit utilizing - = 0.2 is A^ n , The value 
A(£)=pj at £j=l for Geu is set by another duality- 
motivated constraint. R is defined as the ratio of 
deep-inelastic longitudinal and transverse structure 
functions. In the elastic limit: 



Rn(x = 1;Q 2 ) = (4M 2 r /Q 2 ) x (G 2 Er 
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Experiment 


M A 

(published) 


AM A 
new-old 


Miller - D - ANX 8 2,77,73 
Bafcer - D - BiVisi 
Kitagaki - D - FNAL S3 
Kitagaki — D — BNL 9a 


1.00 ±0.05 
1.07 ±0.06 
1 05+ 012 
1 070+ 040 

l.U/U_ 045 


-0.030 
-0.028 
-0.025 
-0.036 



Table 1: Ma GeV/c 2 values published by 
neutrino-deuterium experiments and updated cor- 
rections AMa when re-extracted with updated 
vector form factors. 



Figure 2: The constraint used in fitting G En stip- 
ulates that G 2 En /G 2 Mn = G 2 Ep /G 2 Mp at high £ 

The solid line is ,^, Bp , and |^ Ep I , and the short- 
dashed line is j§ff-\ and ]uff\ ■ 

For inelastic scattering, as Q 2 — > oo, i?„ = 

If we assume quark-hadron duality, the same is for 

the elastic form factors at £=1 ( Q 2 — > oo) 

En/ u Mn ~ u Ep/ Lj Mp 

The new form factors Ge p , Gm p I ' Hp, GMn/Hn, 
and Ge,i are plotted in Figure 1 as ratios to the 
dipole form G^. as ratios to the dipole form fac- 
tor, Gp. An(£) is not needed for Gm p as it is 
very close to 1.0. For Ge p it yields a correction 
of 1% at low Q 2 (because it is required to agree 
with the fits of Arrington and Sick[3] which in- 
clude two photon exchange and Coulomb correc- 
tions. For Geu and Gmu it is used to impose 
quark-hadron duality asymptotic constraints. Fig- 
ure 2 shows plots of -M^r and (for the ^ = 
at £ = 1 case). The long-dashed line is a simple 
quark-hadron duality model [4] . 

Using the updated vector form factors, we per- 
form a complete reanalysis of published neutrino 
quasielastic [6] data on deuterium. (Because of un- 
certain nuclear corrections, neutrino data on heav- 
ier nuclear targets are not used.) We extract new 
values of Ma, and updated values of Fa(Q 2 )- The 
average of the corrected neutrino measurements is 
Ma = 1.0156 ± 0.0278. This is to be compared 
to the average value of 1.014 ± 0.016 extracted 
from pion electroproduction experiments [5] after 
corrections for hadronic effects. The average of the 



two average values is 1.0144 ± 0.0136. The new 
form factors Ge p , Gm p /^ p , Gmu/ 'tin, and Geu 
are plotted in Figure 1 . 

For deep-inelastic scattering, the vector and axial 
parts of the inelastic structure functions W% are 
equal. Local quark-hadron duality at large Q 2 
implies that the axial and vector components of 
-^elastic a j SQ gq Ua j^ yielding: 

rp (n 2, 12 (G V E ) 2 (Q 2 )+T(GUQ 2 )) 2 

[fAKQ )A2=V2\ = (TTt) ' 

(1) 

where G V E {Q 2 ) = G Ep (Q 2 ) - G En (Q 2 ) and 
G V M (Q 2 ) = G Mp (Q 2 ) - G M n{Q 2 ). We do a du- 
ality based fit to the updated values of the axial 
form factor Fa {Q 2 ), including pion electroproduc- 
tion data. Here the fit function is a sum of La- 
Grange polynomials, A a FA , multiplied by G^(Q 2 ) 
(with M A = 1.015). 

F A (Q 2 ) = A FA (0 xG£(Q 2 ). 

We impose the constraint A FA (£i = 0) = p\ = 
1.0. We also constrain the fit by requiring that 
A a FA {0 yield F A (Q 2 ) = F A {Q 2 ) A 2=V2 for £ > 
0.9 (Q 2 > 7.2(GeV/c) 2 ). Figure 3(a) shows val- 
ues of Fa{Q 2 ) extracted from neutrino-deuterium 
experiments divided by G^(Q 2 ), with Ma = 
1.015. Figure 3(b) shows values of Fa(Q 2 ) ex- 
tracted from pion electroproduction experiments 
divided by G^(Q 2 ). These pion electroproduc- 
tion values can be directly compared to the neu- 
trino results because they are multiplied by a fac- 
tor F A (Q 2 ,M A = 1.0U)/Fa{Q 2 ,M a = 1.069) 
to correct for AMa = 0.055 originating from 
hadronic effects [5]. The solid line is our duality 
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Figure 3: (a) Fa(Q 2 ) re-extracted from neutrino-deuterium data divided by Gp(Q 2 ) (with Ma = 1.015). 
(b) Fa{Q 2 ) from pion electroproduction divided by Gp(Q 2 ) (corrected for for hadronic effects[5]). The 
solid line is our duality based fit. The short-dashed line is Fa(Q 2 )a2=V2- The dashed-dot line is a 
prediction from a constituent quark model. 



based fit. The short-dashed line is Fa(Q 2 )a2=V2- 
The long-dashed line is Fa(Q 2 )ai=vi- The 
dashed-dot line is a prediction from a constituent 
quark model[7]. Our new parameterizations of 
vector and axial nucleon form factors use quark- 
hadron duality constraints at high momentum 
transfers and maintain a very good descriptions 
of the form factors at low momentum transfers. 
These parameterizations are useful in modeling 
neutrino interactions (e.g. for neutrino oscillations 
experiments). Our predictions for Geu{Q 2 ) and 
Fa{Q 2 ) at high (Q 2 ) can be tested in upcomng 
electron scattering and neutrino experiments at Jef- 
ferson Laboratory and at Fermilab (MINERvA). 
The final parameters are given in Ref. [1]. 
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